Evidence for Suppression of Superconductivity by Spin Imbalance in Co-Al-Co single 

electron transistors 
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Spin imbalance is predicted to lead to suppression of super- 
conductivity. We report phenomena manifesting this effect 
under spin-polarized quasiparticle currents in ferromagnet- 
superconductor-ferromagnet single electron transistors. The 
measured superconducting gap as a function of magnetic field 
reveals a dramatic decrease when the magnetizations of the 
two leads are in opposite orientations. The effect of suppres- 
sion increases with increasing voltage but decreases at ele- 
vated temperatures. The possible explanations for these de- 
pendences are given. This method may render it applicable to 
control superconductivity at low temperatures and low fields. 

PACS numbers: 74. 25. Ha, 75. 75. -fa, 85.75.-d, 85.35. Gv 



The interplay between superconductivity and mag- 
netism has been a topic of interest for many years. Of 
particular interest are recent experiments on ferromag- 
net/superconductor (FM/SC) junctions with high Tc ||] 
in which a decrease of the supercurrent by nonequilib- 
rium spin density was demonstrated. Furthermore, the- 
oretical studies [HJ^ also indicate that spin imbalance 
in a superconductor can lead to suppression of super- 
conductivity. In a double tunnel junction containing a 
normal metal or a superconductor sandwiched between 
two ferromagnets, both injection of polarized current and 
spin accumulation are possible, and thus provide an ideal 
test ground for the theory. When the magnetic moments 
of the two ferromagnetic leads are in opposite orienta- 
tions (referred to as the antiferromagnetic (AF) align- 
ment), the difference in the number of majority and mi- 
nority spins in the central electrode brings in a chemi- 
cal potential difference. In a normal metal, this poten- 
tial difference is superconductors — PeV/2, where P 
is the polarization of the ferromagnetic leads and V is 
the voltage across the sample. In a superconductor, this 
difference gives rise to pair breaking, same as the Zee- 
man effect does to superconductivity in the paramagnetic 
limit. In this paper, we report the first direct observation 
of the superconducting gap suppression using Co/Al/Co 
(FM/SC/FM) double tunnel junctions. The effect can 
be turned on and off by manipulating the mutual orien- 
tations of the magnetic moment of the two Co leads. We 
further discuss some relevant time scales, including the 
spin relaxation time and the tunneling time, which set 
the criterion for spin accumulation in the superconduc- 



tor. 

The inset of Fig. 1 shows an SEM picture of a measured 
sample and its biasing circuit. The samples were fabri- 
cated by standard electron-beam lithography techniques 
and by the shadow evaporation method. A thin native 
AI2O3 layer between the Al island and the Co electrodes 
acted as a tunnel barrier. Aluminum is a good candi- 
date to serve as the superconductor for this purpose, not 
just for its long spin lifetime Q| allowing the full range of 
spin effects to be studied, but also for its high quality na- 
tive AI2O3 barrier which was shown to have no spin- flip 
tunneling processes |^] . Our electron-gun deposited Al is- 
lands have a superconducting transition temperature of 
about 2K. The samples were measured using a dilution 
refrigerator and the magnetic field was applied along the 
Co leads, i.e. parallel to the long edge of the Al island. 
Because of the small thickness («25nm) of the Al island 
and of the alignment with the magnetic field, the critical 
magnetic field Hc\\ for the Al islands was 21kGs. Exper- 
imentally, we found that the coercivity field for the two 
Co-leads was smaller than 2kGs, well below i?cth ^^"^ 
thus has negligible influence on the superconductivity of 
the Al island. A necessary condition for the imbalance is 
that the spin relaxation time ts is longer than the time 
tt between two successive tunneling events. A rough es- 
timate of Tt is e// 1^, which, as will be shown below, is 
about O.l^I ns. This is shorter than the reported spin 
relaxation time Ts of about 10ns in Al films making 
the spin imbalance possible. 

The current-voltage {IV) characteristics are similar to 
those for a normal metal-superconductor-normal metal 
(N/SC/N) double junction structure. Because the tun- 
nel junction is small, the charging energy Ec associated 
with it is large 0. The current-voltage characteristics 
measured at low temperatures displayed, in addition to 
the superconducting gap 2Ao/e, a pronounced Coulomb 
gap 2Ec/e. This Coulomb gap can be suppressed to a 
minimum by tuning the voltage applied on the gate elec- 
trode, which is coupled to the central island via a capac- 
itor Cg. From a plot of the zero field IV characteristics 
as a function of gate voltages {i.e. the "Coulomb Block- 
ade parallelogram" ||]), we estimated a superconducting 
gap Ao of about 250 ^eV, and a charging energy Ec of 
100~150 /xeV. 

Throughout this study, the samples were symmetri- 
cally current-biased with respect to the ground. We ap- 
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plied a parallel magnetic field to decrease the size of 
the superconducting gap, and monitored this decrease 
by measuring the voltages. As shown in Fig. 1, when the 
field was ramped from +25kGs to -25kGs, we found a 
sharp voltage drop at a field of around -l.SkGs; When 
the field was reversed and ramped up from -25kGs to 
around +1.5kGs, then the voltage dropped again, yield- 
ing a symmetric V{H) pattern with respect to zero field. 
Because of the current-bias circuitry, the change in the 
current was less than 1% regardless of the large voltage 
change. This hysteresis was reproducible, and similar re- 
sults were found for other samples. This behavior is very 
similar to the hysteresis seen in tunneling magnetoresis- 
tance (TMR) of a ferromagnet-insulator-ferromagnet sin- 
gle tunnel junction (see for example Ref. [^). By token 
of this similarity it is reasonable to argue that a field of 
\H\ = 25kGs is sufficiently large to saturate the magne- 
tization of the two leads and to make them parallel to 
the applied fields, whereas at small fields {H — ±1.5kGs 
for the sample of Fig. 1) they are in AF alignment. The 
presence of the two dips in voltage in the AF alignment 
thus demonstrates the suppression of superconductivity 
of the Al island, due to spin imbalance. 

Figure 1 also shows V{H) curves taken at various bias 
currents. By taking the V{H) curves at various bias cur- 
rents one can recover the IV characteristics, as displayed 
in Fig. 2, for three different cases: (I) When the island 
is in the superconducting state with no applied magnetic 
field. (II) When the island is in the normal state and the 
two leads are in the ferromagnetic alignment. (Ill) When 
the island is in the superconducting state and the two 
leads are in the AF alignment. For the case (I), the ob- 
tained IV characteristic can be fitted by a simulation 
of a N/SC/N single electron transistor at gate voltage Vg 
corresponding to the minimum Coulomb blockade thresh- 
old {i.e. Vg — e/2Cg). From this fitting, one can obtain 
a good estimate of the sample parameters, including the 
capacitances and the resistances of the two tunnel junc- 
tions and, most importantly, the zero-field superconduct- 
ing gap Aq of the Al island. The fitted Ao value is about 
260 /ieV, agreeing well with the gap estimate from the 
plot of "Coulomb Blockade parallelogram" . The fitted 
junction parameters can then be checked by a fitting to 
the IV characteristics obtained for the case (II), that is, 
for an all- normal- metal (N/N/N) transistor. The param- 
eters for this particular sample are as follows (with sub- 
scripts S and D stand respectively for source and drain) : 
Rs = i?c=237kW, Cs=360aF, Cr,=300aF, Cg=0.47aF, 
and Vc,=160mV. A temperature of 250mK is used for this 
fitting; this is reasonable considering the heating due to 
eddy current induced by the sweeping field. At a voltage 
interval, the calculated current for the case (II) is smaller 
than the measured ones; this may be attributed to tun- 
neling processes that are not included in the simulation. 
For example, it is known that for Vg = e/2Cg, there is a 
knee at ^ = 2Ec/e = 0.24mV, below which sequential 



tunneling process suppresses the differential conductance 
by a factor of two. However, in reality this feature may 
be smeared probably due to Zener tunneling to the upper 
band. 

With the knowledge of the sample parameters one can 
calculate the IV characteristics atVg — e /2Cg for various 
superconducting gaps Aa. Note that in the AF align- 
ment, the effective drain tunneling resistance Rda in- 
creases by a factor of (1+P^)/(1— P^). For the case of the 
Co electrodes the reported P value is about 0.4. |Q Tak- 
ing into account this correction, we plot the calculated 
curves together with the IV curve for the cases (I)'^(III) 
in Fig. 2. The cross points between the IV curve for the 
case (III) and the calculated curves for various supercon- 
ducting gaps Aa in the AF alignment give a relationship 
between the voltage V and the corresponding Aa, as 
depicted in Fig. 3(a). Theoretically the dependence of 
Aa on V is determined by combining the Aa(<5^) and 
dfj,{eV) functions, and self-consistent calculating the val- 
ues of Sfj,, eV and Aa. |§ The calculated Aa/Ao(V") 
curves are shown in Fig. 3(a). While the theory is based 
on the assumption that majority spins can accumulate 
in the island, this imbalance of spins is a nonequilibrium 
process. In view of this, our analysis, which assumes a 
-Rda using a P value of the ferromagnetic leads, could 
be an oversimplified approach. Nevertheless, this analy- 
sis should give a quantitative picture of the Aa/Ao(V^) 
dependence. 

Our analysis (depicted in Fig. 3(a)) with P=0.4 shows 
that at V < 2Ao/e the ratio Aa/Ao(V^) increases 
with voltage and a maximum value appears at Vpeak ~ 
2.4Ao/e. To see the effect of the P value on the 
Aa/Ao(F) dependence, we further test a P value of 0.6, 
which corresponds to a higher Rda- As shown in Fig. 
3(a), because of this increase in resistance, the Vpoak value 
shifts toward 2Ao/e. Furthermore, for both P values 
at V below Vpcak we found suppression of Aa to about 
Aq{V) and a lowering of Aa with decreasing V. This 
behavior is even pronounced for sample #2. Ideally, at 
low temperatures and low voltages, absence of quasipar- 
ticle tunneling should lead to a Aa/Ao(T^) value close to 
1. Various origins can give rise to excess current, which 
are not taken into account in our simulation. This excess 
current, in our analysis, can give a lower A a value than 
what is expected. One possible source for such excess 
current is the co-tunneling process. In this process, 
the excess current increases at voltages close to Aa/b, 
and has stronger infiuence for IV curves with smaller 
Aa values, and is thus responsible for the decrease of the 
Aa/Ao value at lower voltages. On the other hand, for 
V > Vpcak , the calculated Aa/Aq value is smaller than 
the experimental ones. This discrepancy probably arises 
from a too long energy relaxation time in our small Al 
island, which can smear the effect of spin imbalance. 

We have also measured the V{II) curves (as those in 
Fig. 1) at various temperatures. By using the same anal- 
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ysis technique as above, we traced out separately the 
Ao(T) and Aa(T) dependences, and plotted the ratio 
Aa/Aq as a function of temperature. As shown in Fig. 
3(b), for T < 0.6K this ratio increases with temperature. 
In this temperature range, both the energy relaxation 
time and the spin relaxation time do not change much 
with temperature. Our result, shown in Fig. 3b, may 
be explained by the fast relaxation of domains in the Co 
leads, as we discuss below. 

In the TMR experiments, the two ferromagnetic elec- 
trodes are usually made of two different materials with 
different coercivities. In our case, the AF alignments may 
be ca used by random domain distributions in the two Co 
leads 
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The junction area is about 80nmx80nm, 
which is close to the reported domain size of a thin Co 
wire Therefore, the measured hysteresis shown in 

Fig. 1 is simply a sign of magnetization reversal in prob- 
ably a single domain. The magnetic field was swept at 
a rate of 2.5kGs/min; this is a crucial condition for the 
observed strong suppression of the superconducting gap. 
As shown in Fig. 4(a), a smaller sweep rate results in 
a smaller suppression effect. This is a signature of a 
lesser spin imbalance and indicates imperfect AF align- 
ment of the two leads. We note that this sweep rate de- 
pendence is completely reproducible, and that similar re- 
sults were found for several other samples. The fact that 
the amount of spin imbalance is sweep rate dependent 
indicates that the magnetization reversal is a nonequilib- 
rium process. All magnetic systems that exhibit hystere- 
sis are expected to relax slowly toward the free energy 
minimum. Indeed, when the field is swept from -|-25kGs 
to -l.SkGs (at which Aa is at a minimum) at a rate of 
2.5kG/niin and pause, then the value of Aa decrease with 
time. Moreover, as shown in Fig. 4(b), the rate of this 
decrease is higher at higher temperatures. The relaxation 
time follows a simple Arrhenius dependence on temper- 
ature lll^ , indicating a thermal activated relaxation of 
magnetization in individual domains. The activation en- 
ergy is about 140K. This relaxation can account for the 
lessening suppression of A a at high temperatures. 

In conclusion, the superconducting gap of a small Al 
island, when incorporated in a single electron transistor 
structure in which the two electrodes were made of Co, 
was found to be strongly suppressed by polarized quasi- 
particle injection from the ferromagnetic leads if the two 
leads were in AF alignment. The effect appears to be 
caused by pair breaking associated with spin-imbalance 
in the superconducting island. At elevated temperatures, 
faster relaxation of magnetization of the two ferromag- 
netic leads yields a diminution of spin imbalance. This 
experiment provides a new method of controlling super- 
conductivity. 
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FIG. 1. Measured voltages as a function of applied mag- 
netic field for several selected bias currents: from bottom, 
0.05, 0.20, 0.35, 0.50, 0.60 and 0.75 nA. The inset shows an 
SEM image of one sample. The white object in the center is 
the Al island, and the scale bar at the bottom is lOOnm. 

FIG. 2. Bold curves: the IV characteristics as recovered 
from Fig. 1 using the voltages and the corresponding bias 
currents at B = ±25kGs (leftmost curve), the maximum volt- 
age at B « (rightmost curve), and the voltages at the dips 
at i? « +1.5kGs (middle curve). The two solid curves are 
calculated respectively for the superconducting and normal 
states, and the dashed curves are calculated for various Aa 
values (shown in the figure are 220,200, 180//eV), as described 
in text. The insert shows measured large-scale IV curves at 
B = and B = 25kGs and the calculated IV curves, which 
confirms the obtained sample parameters. 

FIG. 3. (a) Aa/Ao as a function of voltage, as obtained 
from Fig. 2. Calculated curves for T =250mK are also plotted 
for comparison, (b) The ratio Aa(T)/Ao(T) as a function of 
temperature for sample #2. The suppression (and thus the 
imbalance) decreases with temperature due to a decrease in 
lead's magnetization. 

FIG. 4. (a) The ratio Aa/Aq as a function of the magnetic 
field sweep rate. The value of Aa decreases with increasing 
sweep rate and saturates at rates > 2.0kGs/min. (b) The 
relcixation time of lead magnetization as a function of inverse 
temperature which shows a thermally activated relaxation be- 
havior. 
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